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ABSTRACT 


An  axi symmetr i c  numerical  model  of  a  liquid  spray 
is  presented.  The  model  consists  of  coupled  sets  of  partial 
differential  equations  for  the  gas  phase  and  ordinary  dif¬ 
ferential  equations  for  the  liquid  phase.  A  numerical  method 
based  on  the  MAC  method  is  presented  and  calculations  show 
good  agreement  with  experiment.  The  effect  of  a  ceiling  on 
the  entrainment  performance  of  a  spray  is  found  to  be  negli¬ 
gible.  The  examination  of  the  calculated  pressure  field 
shows  that  pressure  does  not  play  an  important  role  in  the 
flow  and  may  be  ignored.  This  is  in  agreement  with  experiment 
and  contributes  additional  validation  to  existing  models 
which  ignore  this  effect. 
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1.  INTRODUCTION 


The  application  of  liquid  sprays  is  made  in  many 
different  fields.  Examples  of  a  few  are  fire  suppression, 
combustion,  ventilation, and  the  dispersion  of  heavy,  possibly 
toxic  or  flammable,  gases-  Interest  at  V K I  in  the  applica¬ 
bility  of  water  sprays  for  use  as  water  curtains  for  the  dis¬ 
persion  of  heavy  gases  has  mandated  a  need  for  better  predic¬ 
tion  methods  for  the  analysis/design  of  water  sprays. 

The  conf i guration  for  the  spray  to  be  analysed  in 
this  study  is  shown  in  figure  1.  Liquid  particles  are  ejected 
through  a  nozzle  into  an  environment  which  is  at  rest  in  the 
absence  of  the  spray.  The  aerodynamic  drag  acting  on  the 
liquid  particles  results  in  a  loss  in  momentum  of  the  particles. 
Since  momentum  must  be  conserved,  the  momentum  loss  is  trans¬ 
lated  into  a  momentum  gain  by  the  fluid,  causing  the  fluid  tu 
be  set  in  motion.  Air  is  thereby  entrained  into  the  spray 
(in  much  the  same  way  as  a  jet)  resulting  in  a  two  phase  plume. 

Modeling  of  this  phenomenon  has  been  done  in  the 
past.  A  one  dimensional  model  (Refs.  3,5)  has  been  developed 
at  VKI  and  has  shown  good  results  though  its  appl i caoi 1 i ty 
for  the  inclusion  of  boundary  effects  is  limited.  Two  dimen¬ 
sional  models,  existing  in  the  literature  (Refs.  1,4,7),  are 
much  more  able  to  handle  the  effects  of  boundaries.  In  this 
light,  it  is  found  desirable  to  have  such  a  two  dimensional 
model  at  VKI  to  complement  the  range  of  applicability  of  the 
existing  one  dimensional  model.  Thus,  this  project  concerns 
the  implementation  of  an  axisymmetric  numerical  model  into  a 
computer  code.  The  basic  approach  is  similar  to  that  followed 
in  references  1  and  4. 

In  this  report  the  model  used  is  outlined,  the 
numerical  method  for  the  solution  of  the  model  equations  is 
described,  and,  finally,  the  results  presented  and  discussed. 
The  implemented  computer  coua  is  included  as  an  appendix  to 
this  report. 
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2.  THE  SPRAY  MODEL 

2.1  Introduction 


The  model  of  the  spray  consists  of  two  distinct 
sets  of  equations,  one  set  governing  the  gas  phase  and  another 
governing  the  liquid  phase.  Linkage  between  these  two  sets  of 
equations  accounts  for  the  following  two  physical  phenomena  : 
the  first  concerns  the  momentum  transfer  between  liquid  phase 
and  gas  phase,  while  the  second  involves  the  modification  of 
the  particle  trajectories  by  the  motion  of  the  gas  phase. 

The  gas  phase  may  be  modelled  as  a  continuum  using 
an  Eulerian  approach  while  the  liquid  phase,  using  a  uagrangian 
approach  is  modelled  by  considering  a  finite  set  of  particles 
of  varying  size/initial  trajectory. 

2 . 2  The  gas  phase  model 

2.2.1  Gas^  £h£se  equate  on^ 

Since  the  gas  phase  occupies  the  most  significant 
portion  of  the  flow  domain,  it  may  be  treated  as  a  continuum. 
Making  the  standard  assumptions  of  i ncompressi bi 1 i ty  and 
Newtonian  fluid  the  following  form  of  the  Navier  Stokes 
equations  may  be  used.  Here,  due  to  the  axisymmetric  nature 
of  the  problem,  the  equations  are  written  in  cylindrical 
coordinates  : 
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where  Vr,  Vz,  P  are  the  radial  and  axial  component  of  velocity, 
and  the  pressure  respecti vely ;  R  and  z  are  the  radial  and  axial 
coordinate  and  p  and  v  represent  the  fluid  density  and  kine¬ 
matic  viscosity.  The  terms  Fr  and  Fz  are  the  volume  force 
terms  which  account  for  the  contribution  of  momentum  per  unit 
volume  from  the  liquid  phase. 


2.2.2  Bounda_ry_cond  i  ti  onj>  for  j:he  jjas^  phjase 

Considered  here  are  three  types  of  boundaries;  a 
single  axis  of  symmetry  boundary,  free  boundaries  and  wall 
boundaries.  These  are  shown  in  figure  1. 

Velocity  boundary  conditions 
Axis  of  symmetry 

At  the  axis  of  symmetry  the  radial  component  of 
velocity  must  be  equal  to  zero  because  of  mass  conservation 
and  the  symmetry  assumption.  Symmetry  also  requires  that  there 
be  no  shear  stress  in  the  axial  direction,  requiring  the  deriva¬ 
tive  of  the  axial  velocity  with  respect  to  the  radial  coordinate 
to  be  zero. 
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Vr  =  0  (2.4) 


dR 

Mall  boundary 

Along  the  wall  boundaries  the  usual  conditions  for 
an  impermeable,  non-slip  wall  are  used  : 

Vr  =  V2  =  0  (2.6) 

Free  boundary 


In  the  flow  domain  considered  the  free  boundary  is 
taken  to  be  far  from  the  spray  domain.  With  this  assumption 
the  following  approximate  conditions  may  be  applied.  The 
velocity  tangential  to  the  boundary  is  small,  and  may  be  assumed 
equal  to  zero.  Applying  continuity,  the  boundary  condition  for 
the  component  of  velocity  normal  to  the  boundary  may  be  obtained. 

Vertical  free  boundary 
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Horizontal  free  boundary 
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Pressure  boundary  conditions 

The  pressure  boundary  conditions  may  be  derived 
from  the  momentum  equations  2.2  and  2.3.  Since  most  gases 
have  small  viscosity  it  is  permissible  here  to  ignore  the 
viscous  terms  of  the  momentum  equations  in  deriving  the  pres¬ 
sure  boundary  conditions. 

The  resulting  derivation  shows  that  for  the  wall  and 
axis  of  symmetry  boundaries,  the  normal  derivative  of  the 
pressure  is  equal  to  zero. 

The  resulting  boundary  condition  on  the  free  boundary 
is  a  bit  more  complicated  using  this  approach.  However,  ma k i ng 
use  of  experimental  observation  which  has  shown  that  variation 
of  the  pressure  field  is  small  throughout  the  flow  domain,  the 
pressure  at  the  free  boundary  is  assumed  to  be  constant. 

Wall  or  axis  of  symmetry 


a  N 

Free  boundary 


(2-11) 


P  =  P 


ref 


(2.12) 


2 . 3  The  particle  phase  model 
2.3.1  f.ailtlcle_e£ua  t^ons 

The  liquid  phase  is  modelled  as  a  distribution  of 
droplets  of  distinct  sizes  and  trajectories.  Using  a  lagrangian 
approach  individual  particles  are  followed  from  injection  until 
hitting  the  floor.  By  considering  the  number  of  particles  of 
similar  sizes/trajectories,  the  magnitude  of  droplet-gas 
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momentum  exchange  throughout  the  flow  field  may  be  determined. 
This  approach  assumes  that  the  particles  do  not  interact  with 
each  other,  i.e.,  no  collisions  or  particle  break  up.  Making 
the  additional  assumptions  that  the  particles  are  spherical  and 
non  -evaporati ng ,  the  following  equations  of  motion  may  be 
wr i tten . 


(2-13) 


draft) 

dt 
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dzo(t) 

dt 
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(2.14) 

(2.15) 

(2-16) 


where  ur>  uz  are  the  particle  radial  and  axial  velocities  and 
ro  and  zo  are  the  radial  and  axial  position  of  the  particles. 
"m“  is  the  mass  of  the  particle  calculated  as  follows  : 


(2.17) 


where  is  the  density  of  the  liquid  and  D  is  the  particle 
diameter.  f  r  and  f  are  the  radial  and  axial  component  of 
drag  force,  related  to  Fr  and  Fz  in  equations  2.2  and  2.3  and 
are  calculated  as  follows  : 


nDv  p 

fr  =  CD  Re  __  (ur-Vr) 
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f*  ■  co  — -  <Vvz> 


(2.13) 
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(2.19) 


Re  is  the  Reynolds  number  defined  as 


Re  = 


y(yV2)2+(Vvr)2  o 

V 


(2.20) 


The  drag  coefficient  is  calculated  using  a 
standard  form  fit  for  the  drag  coefficient  of  a  sphere  : 

Cn  =  —  +  — +  .4  (2.21) 

u  Re  1+y  Re 

2.3.2  Ini_t2al_  condiyons  f  or_t he_pa/t i  cl  e  equations 

The  initial  conditions  are  derived  from  the  proper¬ 
ties  of  the  spray  nozzle.  By  definition  the  coordinate  system 
is  fixed  at  the  nozzle  so  the  initial  conditions  for  equations 
2.15  and  2.16  are  simply  : 

r0(0)  =  z0(0)  =  0  (2.22) 

The  initial  conditions  for  the  force  equations  are 
derived  by  considering  the  particle  ejection  velocity  from  the 
nozzle  and  the  initial  angle  of  ejection  (see  figure  2)  : 

ur(0)  =  U0sine 

uz(0)  =  U0coso  (2.23) 


Ug  is  calculated  from  the  volume  flow  of  the  nozzle,  )  : 


U 


o 


(2.24) 
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where  do  is  the  inside  diameter  of  the  nozzle.  The  initial 

angle  of  ejection  of  the  particle  must  be  between  0  and  e 

0  and  dft  are  available  from  the  manufacturers  data, 
max  u 


2.3.3  £a.rtic]_e_s2ze 

Particle  size  i'  calculated  by  using  the  following 


expression 


D 


2/3 


(2.25) 


where  is  a  nozzle  design  parameter  and  aP^  is  the  water 
delivery  pressure.  These  values  may  be  determined  from  experi¬ 
ment  or  from  manuf acturers  data. 


2 . 4  Nondimensio nalization  of  the  model  equations 

To  allow  the  resulting  program  to  be  used  easily 
with  any  system  of  units,  the  model  equations  are  nondimension- 
al i zed . 


Velocities  are  nondimensionalized  by  the  particle 
ejection  velocity 


A 


V 


r 


(2.26) 


Pressure  is  non 
and  the  particle  ejection 


dimensional i zed  by  the  gas  density, 
velocity  squared 


(2.27) 


Lengths  are  nond i me  ns i o na 1 i zed  by  the  spray  height,  Hf,  (see 
figure  1) 


A  o  A  7  A  ro  A  zo 

R  =  — ,  z  =  — ,  r0  =  — ,  z0  =  — 

H*  Ht  Ht  Ht 


(2.28) 


Time  is  nondi mens i o na 1 i zed  by  the  spray  height  divided  by  the 
ejection  velocity 


A 

t 


(2.29  ) 


The  momentum  source  terms  are  nond i me  ns i ona 1 i zed  by 
the  ejection  velocity  squared  divided  by  the  spray  height 


A 

r 


(2.30) 


Thus  in  the  momentum  equations  (2.2)  and  (2.3),  the 
viscosity  is  replaced  by  the  reciprocal  of  the  "Reynolds  number 
of  no nd i mens i ona 1 i za t i on " 


JL  = 

RE  U0H^ 


(2.31) 


The  nondi me nsionalization  of  the  other  model  equations 
is  straightforward.  From  this  point  on,  th  ■  model  equations  are 
assumed  nondi, r.ensional  ,  and  the  hat  (a)  is  neglected. 
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3.  THE  NUMERICAL  METHOD 

3 . 1  Introduction 

It  is  difficult  indeed  to  conceptualize  a  numerical 
method  for  the  simultaneous  solution  of  both  sets  of  equations. 

A  more  natural  method  of  solution  involves  iterative  solution 
of  these  equations.  Such  an  iterative  solution  is  outlined 
by  the  flow  chart  in  figure  3. 

The  specific  numerical  method  to  solve  each  set  of 
equations  differs  due  to  their  nature  (ODE,  PDE)  and  are 
considered  in  appropriate  sections. 

3 . 2  Numerical  solution  of  the  particle  equations 

The  solution  of  the  set  of  ordinary  differential 
equations  (ODE)  governing  the  particle  phase  is  made  by  a 
s imu 1 taneous  fo urth  order  Runge-Kutta  technique  which  is  well 
documented  in  the  literature  (Ref.  9). 

3 . 3  Numerical  solution  of  the  gas  equations 

The  numerical  scheme  for  the  solution  of  the  gas  phase 
equations  is  based  on  the  MAC  method  (Refs.  8,10,11). 

The  gas  phase  equations  are  recast  into  the 
following  form 
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These  equations  are  written  in  unsteady  form  as  the 
MAC  solution  procedure  is  an  iterative  one  in  which  the  steady 
state  solution  is  the  desired  solution.  To  be  noted  here  is 
the  exchange  of  the  continuity  equation  for  a  Poisson  equation 
for  the  pressure.  Contained  in  this  Poisson  equation  is  the 
variable  D  which  is  a  dilation  term  representing  the  amount  of 
continuity  existing. 


0 


l  3(VrR)  3 Vz 

R  J  R  3Z 


(3.4) 


This  term  is  used  to  reduce  the  nonlinear  instabilities  existing 
in  the  numerical  solution  of  the  Na v i er- S to kes  equations 
(Ref.  11).  Continuity  is  solved  implicitly  using  this  term. 
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3.3.1  Discretization  of  the  MAC  method  equations 

An  appropriate  discretization  of  the  MAC  method 
equation  on  a  staggered  grid  (see  Fig.  4)  is  given  below  : 
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Here  the  tilde  ( ~ )  indicates  the  updated  quantities 
(N+l  i teration ) . 

As  can  be  seen,  the  MAC  method  is  a  two  level  scheme 
involving  explicit  solution  of  equations  3.5  and  3.6  for  the 
velocities  of  the  N+l  iteration  and  an  iterative  solution  of 
equations  3.7  for  the  updated  pressure  field. 


Note  the  discretization  of  the 


term. 


is  set 


equal  to  zero  in  an  attempt  to  force  continuity  to  exist  at  the 
N+l  time  step.  At  steady  state  the  time  derivative  term  disap¬ 
pears  and  continuity  will  be  satisfied,  thereby  solving  the 
original  problem. 
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This  method  was  applied  to  the  solution  of  laminar 
flow  in  a  pipe  as  a  test  case.  Solutions  were  obtained  for 
low  Reynolds  numbers  (Re  <  100)  but  instabilities  were  observed 
for  higher  Reynolds  numbers. 

Because  the  solution  of  the  gas  equations  involves 
a  low  gas  viscosity  (v  =  1.5  *  10~5  m2/s)  and  therefore  high 
Reynolds  numbers,  the  numerical  method  has  to  be  modified. 

3.3.2  Mod i f i cati on_to  the  MAC  method 

To  remove  the  instabilities  upwind  differencing  was 
employed.  This  is  not  so  straightforward  for  the  staggered 
grid  used  and  involves  some  additional  approximations.  The 
resulting  discretization  of  the  momentum  equations  is  made 
below.  Here  only  the  convective  terms  are  affected  and  there¬ 
fore  are  the  only  ones  considered. 
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The  pressure  equation  remains  the  same. 
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This  irodi  ri  cation  allowed  a  stable  solution  to  be 
obtained,  though  the  scheme  did  not  allow  the  lower  wall  to 
be  sensed,  and  no  rec i rcu 1  a t i on  occurred  (see  Fig.  5).  There 
another  modification  was  employed.  This  involved  changing 
the  unsteady  term  in  the  momentum  equations  along  the  lower 
and  axis  of  symmetry  boundaries. 
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on  lower  boundary. 
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This  finally  yielded  plausible  results. 

3.3.3  Sol^ut  i  on_procedur  e_f  or_tjie_MAC_meth^od 

The  solution  procedure  for  the  MAC  method,  as 
mentioned  previously,  involves  an  explicit  solution  of  equa¬ 
tions  (3.5)  and  (3.6)  with  the  modifications  for  the  velocity 
field  at  the  N  +  l  iteration  and  an  iterative  solution  of  (3.7) 
for  the  pressure  field  at  the  N+l  iteration.  This  iterative 
solution  was  made  using  a  point  by  point  SOR  method.  The 
overrelaxation  factor  used  was  w  =  1.5,  though  no  optimization 
has  been  done  on  this  parameter. 


la 


3 . 4  Calculation  of  the  source  terms 
in  the  Navier-Stokes  equations 


For  the  calculation  of  the  source  term,  the  particle 
trajectories  are  superimposed  over  the  mesh  used  for  the  solu¬ 
tion  of  the  gas  equations,  (see  figure  6).  To  calculate  the 
source  terms  the  arc  length  of  each  particle  trajectory  in  a 
given  cell  is  determined.  Because  the  velocity  along  this 
trajectory  is  known,  the  time  that  a  particle  spends  in  this 
cell  can  be  determined  : 
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(3.12) 


From  this  time  the  average  number  of  particles  in 
the  cell  is  determined  by  considering  the  fraction  of  the  mass 
flow  through  the  nozzle  assigned  to  this  particular  trajectory, 
fp.  The  average  force  along  the  arc  length  is  determined  and 
then  the  contribution  to  the  momentum  source  term  along  this 
trajectory  in  this  particular  cell  is  found  as  follows. 


Number  of  particles  on  given  trajectory  in  control 

volume  : 
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ra  v 


Vol ume  of  cell 


(3.14) 


simi  1  ar ly  for  Fz  . 

It  should  be  noted  that  Fr  and  Fz  are  not  calculated  in  the  same 
cell  because  of  the  staggered  grid  (see  figure  4). 
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4.  RESULTS 


Results  were  obtained  for  the  consideration  of  the 
effect  of  increase  of  mass  flow  through  the  nozzle  and  the 
effect  of  a  ceiling  on  entrainment  properties  of  the  spray, 
and  finally,  to  examine  the  pressure  field.  The  results  were 
obtained  by  using  the  properties  of  a  Lechler  SZ1  spray  nozzle 
as  input  to  the  program  (see  Table  1). 

4 . 1  Comparison  with  experiment 

First  to  validate  the  program  a  comparison  of  the 
entrainment  properties  of  the  calculated  spray  was  made  with 
the  experimental  results  of  an  unconfined  spray.  These  results 
are  shown  in  figure  7.  Here  the  entrainment  efficiency  is  plot¬ 
ted  versus  the  inverse  of  the  non  dimensional  spray  envelope 
diameter  squared.  The  entrainment  efficiency  is  defined  as  the 
volume  flow  of  air  entrained  into  the  spray  divided  by  the 
volume  flow  of  water  through  the  nozzle.  The  amount  of  air 
entrained  into  the  spray  is  defined  as  follows  : 


Q,  =  2i 


t/2 


Vz(r) 


d  r 


(4.1) 


Results  were  obtained  by  considering  a  spray  nozzle  mounted  on 
a  ceiling  two  meters  above  the  floor,  for  all  three  mass  ; ' ows 
considered.  Although  this  is  not  an  unconfined  spray,  previous 
results  have  shown  little  effect  of  the  presence  of  the  ceiling 
(Ref.  6).  The  axial  positions  where  these  results  were  obtained 
ranged  from  20  percent  of  the  spray  height  from  the  nozzle  to  55 
percent  of  this  distance.  It  was  impossible  to  calculate  these 
values  any  closer  to  the  nozzle  because  of  the  lack  of  resolution 
with  21  x  21  pressure  nodes.  At  axial  positions  greater  than 
this  the  floor  seems  to  have  a  significant  effect  on  the  entrain¬ 
ment,  so  no  comparison  is  made  in  this  area. 
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The  results  displayed  in  figure  7  show  good  agreement 
with  the  McQuaid  correlation  (solid  line)  and  the  i  201  scatter 
of  experimental  data  about  this  line.  This  gives  confidence 
that  the  bulk  of  the  flow  field  is  well  predicted. 

4 . 2  Effect  of  increasing  mass  flow  through  nozzle 

The  effect  of  increasing  mass  flow  is  displayed  in 
figures  8,  9  and  10.  Here  the  velocity  vectors  are  non  dimen- 
sionalized  by  the  ejection  velocity  of  the  particle.  For  the 
determination  of  the  entrained  air  flow  equation  4.1  is  rearranged 

(t/2  V  (r)  <t/2  V  (r) 

Q  =  2tt  U0  1  r  — -  dr  a  2tt  Q  1  r  — -  dr  (4.2) 

a  J  U„  w  J  ll„ 


As  can  be  seen  from  the  results,  the  flow  field  non  dimensio- 
nalized  by  the  ejection  velocity  does  not  change  significantly. 
Thus  the  integral  in  4.2  is  of  the  same  order  of  magnitude  for 
all  three  mass  flows  considered.  It  seems  that  the  air  entrain¬ 
ment  in  the  spray  is  increased  strongly  by  increasing  the  mass 
flow  through  the  nozzle.  This  is  a  physically  observable  result. 


Finally  the  efficiency  of  the  spray  is  examined. 
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Close  examination  of  figures  8,  9  and  10  shows  a  very  slight 
increase  in  the  efficiency  of  the  spray  with  increased  mass 
flow,  reflected  by  a  slight  increase  in  the  above  integral. 

This  result  is  found  to  be  true  for  some  cases,  though  the  in¬ 
verse  can  also  be  true.  This  is  because  higher  mass  flow 
results  in  small  droplets  which  exchange  momentum  with  air  more 
rapidly  resulting  in  higher  entrainment,  but  because  of  the 
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small  droplet  size  the  spray  envelope  contracts  quickly  slowing 
entrainment.  Thus,  for  small  spray  lengths  an  increase  in 
efficiency  is  observed  relative  to  larger  droplet  sizes,  but 
this  is  reversed  as  the  spray  length  is  increased. 

4 . 3  Effect  of  ceiling  on  entrainment 

Results  were  obtained  for  the  case  where  the  ceiling 
is  removed.  These  results  are  shown  in  figures  11  and  12. 
Though  not  shown  a  free  boundary  was  placed  at  a  height  of  two 
meters  above  the  spray  nozzle.  Comparison  with  the  previous 
results  of  corresponding  mass  flow  with  ceiling  included  shows 
little  change  in  the  entrainment  properties  with  or  without  a 
ceiling.  This  result  is  supported  by  experimental  results 
(Ref.  3). 


4 . 4  Examination  of  the  pressure  field 

Along  with  the  velocity  field  each  calculation  also 
gives  the  corresponding  pressure  field.  The  isobars  of  a 
typical  pressure  field  are  shown  in  figure  13.  This  figure 
validates  the  previous  assumption  in  the  derivation  of  the 
boundary  condition  that  at  the  free  boundary  the  pressure  field 
varies  little  and  the  pressure  can  be  set  to  a  reference  value. 

The  area  of  steepest  pressure  gradients  is  near  the 
floor  inside  the  spray  envelope.  However,  the  pressure  gradient 
in  this  region  is  small  compared  to  the  momentum  source  term. 
The  other  area  of  high  pressure  gradient  seems  to  be  at  the 
nozzle.  In  this  region,  again  the  pressure  gradient  is  small 
compared  to  momentum  source  term  inside  the  spray  envelope, 
though  outside  the  spray  envelope  the  pressure  gradient  may 
be  significant.  Better  resolution  is  necessary  to  draw  a  firm 
conclusion. 
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It  can  be  concluded  that  for  the  bulk  of  the  flow 
field  ignoring  the  pressure  gradient  may  be  a  good  approxima¬ 
tion.  This  is  a  particularly  important  result  since  the  pres¬ 
sure  calculations  are  the  most  time  consuming  and  require  the 
most  number  of  statements  in  program.  In  the  future  the 
program  may  be  modified  to  include  only  a  closed  form  expres¬ 
sion  for  a  pressure- 1 i ke  variable  which  forces  continuity  to 
be  satisfied.  This  should  result  in  large  savings  in  computer 
time  and  allow  better  resolution  in  the  mesh. 


5.  CONCLUSIONS 


It  has  been  shown  that  the  axisymmetric  spray  model 
implemented  gives  good  results  for  the  entrainment  of  air  in¬ 
side  the  spray  envelope  for  the  portion  at  the  spray  envelope 
which  contains  most  of  the  entrained  air.  The  bulk  of  the 
surrounding  flow  field  is  thus  considered  to  be  more  or  less 
correctly  predicted. 

The  physically  observable  result  that  increasing  the 
outflow  from  the  nozzle  increases  entrainment  of  air  is  also 
found . 


The  effect  of  the  ceiling  on  entrainment  is  seen  to 
be  negligible. 

Finally,  the  examination  of  the  pressure  field  indi¬ 
cates  that  neglecting  the  pressure  gradient  is  a  good  approxi¬ 
mation  for  the  bulk  of  the  flow  field.  This  result  leads  to 
the  recommendation  that  work  continues  on  the  model  to  incor¬ 
porate  a  closed  form  expression  for  a  pressure-like  variable 
to  force  continuity  to  hold,  while  eliminating  the  need  for  the 
Poisson  equation  for  the  pressure.  This  will  allow  better 
results  since  the  mesh  may  be  made  smaller  due  to  the  saving 
in  computer  time. 

Application  of  this  model  to  a  spray  in  the  upward 
facing  mode  is  straightforward. 
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